Power law Kohn anomaly in undoped graphene induced by Coulomb interactions 
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Phonon dispersions generically display non-analytic points, known as Kohn anomalies, due to 
electron-phonon interactions. We analyze this phenomenon for a zone boundary phonon in undoped 
graphene. When electron-electron interactions with coupling constant /? are taken into account, one 
observes behavior demonstrating that the electrons are in a critical phase: the phonon dispersion 
and lifetime develop power law behavior with /3 dependent exponents. The observation of this 
signature would allow experimental access to the critical properties of the electron state, and would 
provide a measure of its proximity to an excitonic insulating phase. 



I. INTRODUCTION 

The study of electron-electron interactions in 
graphene, a monolayer of carboni, remains one of the 
most interesting open problems in the field and is 
currently a very active area of research^. The Coulomb 
interaction plays a very particular role in this system be- 
cause its low-energy electronic excitations are described 
by massless Dirac fermions. For undoped graphene, 
the vanishing of the density of states at the Fermi level 
implies that the interaction remains truly long-ranged, 
decaying in real space as 1/r, and this is predicted 
to lead to a number of exotic interaction-induced 
phenomena, such as logarithmic rcnormalization of 
many physical obscrvablcs at weak coupling^—, and 
instabilities towards different symmetry breaking states, 
like the excitonic insulator^, at strong coupling. The 
strength of the Coulomb interaction in graphene can 
be characterized by its bare fine structure constant 
j3 = e 2 /evF, where vp is the Fermi velocity and e the 
dielectric constant. A naive estimate yields (3 « 2 for 
suspended graphene, while lower values arc obtained 
in the presence of a substrate. This suggests that 
Coulomb interactions can be relatively important, but 
experimentally their strength is still debated^. 

One of the most striking predicted effects of the 
Coulomb interaction in graphene is that, because the 
Hamiltonian that describes it is scale invariant, some of 
its correlation functions behave like power laws with in- 
teraction dependendent exponents, so that the system is 
effectively in a critical phased. This has also been shown 
recently by renormalization group argument o 10 ' 11 . Sim- 
ilar power laws are also found in Dirac fermion models 
with interactions mediated by effective gauge field o 12 ' 13 . 
Unfortunately most of the usual experimental probes do 
not couple to the correlation functions that display criti- 
cal behavior, making their experimental observation chal- 
lenging. 

In this work we show that a signature of this criticality 
may be accessed experimentally through the dispersion 
relation of a zone boundary phonon, the A\ phonon at 
the K point. The dispersion of this phonon is produced 
mainly by its interaction with the Dirac electrons. With- 
out electron-electron interactions it shows a square root 



cusp at q = qx = ojk/vf (^k the phonon frequency 
at the K point) that crosses over to linear dispersion for 
q >> qx, a feature known as a Kohn anomaly. In this 
work we demonstrate that when the Coulomb interaction 
is included, the phonon dispersion is modified strongly: 
around qx it becomes a power law cusp with exponent 
r]{0), and for q >> qx it crosses over to another power 
law with exponent 770 (/3) . The observation of this strong 
modification of the Kohn anomaly, in principle feasible 
with current experimental techniques^—, would provide 
dramatic evidence of the critical Coulomb interactions in 
this system, and could potentially be used as a much 
needed measurement of their strength j3. This remark- 
able power law Kohn anomaly is similar to the one found 
in some one dimensional systems^. 

The presence of these powers laws can be understood 
in simple terms, while their detailed behavior requires 
an elaborate calculation discussed below. Consider the 
usual low-energy Hamiltonian for graphene around the 
K and K' points 



H = iv F / d 2 r^{a x d x + a y dy)if;, 



(1) 



with a = (r z <7 x , a y ), where the a and r matrices act 
on the sublattice and valley degrees of freedom, respec- 
tively (spin will be accounted for when necessary). The 
chemical potential is set to zero. This Hamiltonian has 
an SU(2) valley symmetry generated by the matrices 
T n = (Tx<r v , Tyijy, t z ), in the sense that the SU(2) ro- 
tation tp — > e' lTn0rl tp leaves the Hamiltonian invariant. 
When the Coulomb interaction 
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(2) 



is included and for /3 greater than some critical value j3 c , 
this system has an instability to an ordered state known 
as the excitonic insulator— where charge imbalance be- 
tween sublattices, i.e. an expectation value of ip'cr z ip, de- 
velops. This instability is reflected in the corresponding 
susceptibility a z ipip^ a zip') at ui = 0, which develops 
a power law q 110 with an interaction dependent exponent 
that goes to zero for j3 — > j3 c , signaling the onset of the 
excitonic phase^^. Therefore, power law behavior in 
this correlator can be thought of as the weak coupling 



2 



counterpart of the excitonic instability, and experimen- 
tal access to the exponent would allow one to probe how 
close the system is to it. However, this particular suscep- 
tibility is difficult to measure, as it requires a probe that 
couples differently to the two sublattices. 

The charge density wave instability is however not the 
only one that Coulomb interactions can induce. The 
Hamiltonian ([I} admits two other time-reversal invari- 
ant masses t x <j x and t v <j x , and an instability that de- 
velops an expectation value for either of them may pro- 
ceed in the same way. These order parameters corre- 
spond to a bond density wave order known as the Kekule 
distortion^. It can be shown that the Kekule and CDW 
masses, M n = (r x a x ,T y a x ,<7 z ) transform like a spin 1/2 
under the valley symmetry, and since the Coulomb inter- 
action does not break this symmetry, the three instabili- 
ties are in fact equivalent: they have the same weak cou- 
pling power law susceptibility with the same exponents. 
Since, as we will see, the A\ electron-phonon vertex cor- 
responds to the Kekule mass, the phonon self-energy is 
proportional to the Kekule susceptibility. Therefore we 
expect power law behavior in the phonon dispersion and 
lifetime. To see this, however, the computation of the 
full uj dependent susceptibility is required. In the re- 
mainder of this paper we discuss the electron-phonon 
coupling in graphene and the computation of the phonon 
self-energies with the aim of establishing precisely where 
the signatures of critical behavior are to be found. 



II. PHONONS AND KOHN ANOMALIES 

The phonon spectrum of the honeycomb lattice con- 
sists of six phonon branches, four in-plane and two out- 
of-planc. Each of these phonons may couple to elec- 
trons near either Dirac point if it has momentum close 
to zero (a T point or zone center phonon), which scat- 
ters electrons within each valley, or if it has momen- 
tum close to K or K' points (a zone boundary phonon), 
in which case it produces intervalley scattering. The 
strength of the electron-phonon coupling (EPC), how- 
ever, depends on how the particular displacement pattern 
of that phonon modifies the hopping integrals between 
atoms. Two modes have displacements that produce a 
significant EPC. The first of these is the phonon branch 
of highest energy at the T point, the E 2 phonon. The 
second is the Ai branch at the K and K' points (also the 
highest branch) . This is a lattice distortion with a super- 
cell of six atoms, whose displacement pattern is obtained 
by taking linear combinations of the displacements at K 
and K', and is shown in the inset of Fig 2J These two 
combinations couple to electrons exactly in the same way 
as the two components of the Kekule distortion 



H, 



e—ph,K 



= F h 



d 2 r^ {Miukx + M 2 u K2 )i), (3) 



with Fk = 3dt/da. For this reason this phonon is also 
known as the Kekule phonon^. The fact that the E 2 



and A\ phonons are the most predominant is confirmed 
by Raman spectroscopy in pristine graphene, where two 
main peaks are observed 2 ^; the G peak corresponds to 
E 2 phonons, while the 2D peak is a second order process 
involving two A\ phonons. The Hamiltonian of the A\ 
phonon may be expressed as 



H = J2 



(27T) 



;U} K b{ n h 



(4) 



with creation and destruction operators defined by 



An 



d 2 q 

4w K M J (2tt) 2 



(b 



b L e ~ 



% (5) 



where i = K1,K2, ujk ~ 0.17 eV, A c is the unit cell 
area. For the range of momenta q < 0.25A -1 where 
the Dirac fermion model is applicable^, the dispersion 
of the phonon can be neglected. Indeed phonon band- 
structure computations excluding the effect of electron- 
phonon coupling show a practically fiat dispcrsio: 



.21.22 



this range. A dimcnsionlcss EPC can be defined as 



(6) 



which is estimated to be in the range Xk ~ 0.03— 0. 1 23 ' 24 . 

Due to electron-phonon interactions, phonon disper- 
sion relations are known to develop non-analytic points, 
known as Kohn anomalies, at the largest momenta for 
which the generation of an electron-hole excitation is 
kinematically allowed. This renormalization of the dis- 
persion, as well as the phonon lifetime, can be ob- 
tained from the phonon self-energy S, which enters in 
the phonon Green's function as 



G ph {oj,q) = — j- 



2uj k 



2uj K ^{u>,q)' 



(7) 



As anticipated, this self-energy is related to the mass 
susceptibility, which is defined as 



n„ m = <V' t M„v ^M m ip) 



(8) 



because of the form of the coupling given in eq. ©. The 
explicit relation follows from the previous definitions and 
reads 



s = ^n n = ^n 22 . 
2 2 



(9) 



In the absence of electron-electron interactions, the self- 
energy can be computed analytically, and it is given 
by2^5 



^ 1 q) = ^f(vW-^) 1/2 - 



(10) 



Solving for the pole in Eq. ([7]) for small A^, we see the 
dispersion relation is corrected to 

Lu(q) a LJ K + X K /A{v 2 F q 2 - UJ 2 ,) 1 ' 2 (11) 
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which has a square root singularity at qK for q > qx- For 
q < qK the self-energy is purely imaginary, and a finite 
lifetime is obtained. The Kohn anomaly is conventionally 
associated with a linear cusp in the dispersion, which is 
obtained only assymptotically for q >> qx] the full dy- 
namical self-energy should be used in general. Note that 
qK is approximately 2% of the V — K distance in the 
Brillouin zone. The necessity of employing the dynami- 
cal self-energy has been emphasized before 2 ^— , in par- 
ticular in the doped case where the static approximation 
produces poor agreement with experiments 2 ^. 



III. POWER LAW MASS SUSCEPTIBILITY 
AND PHONON DISPERSION 

We will now proceed to compute the general uj 
and q dependent mass susceptibility Tl nm including the 
Coulomb interaction. We will see that it acquires /3 de- 
pendent power law behaviour, a feature that is thus in- 
herited by the A\ phonon. We will employ a ladder sum- 
mation, as it is the simplest approximation that will cap- 
ture any non-analytic behaviour. The ladder summation 
is represented diagramatically in Fig. [TJ Denoting three- 



® 



(a) 





(b) 



FIG. 1. (a) Diagrammatic equation for the three point ver- 
tex (shaded triangle) in the ladder approximation. The cross 
denotes a mass vertex, (b) Response function diagram. 



momenta q = (qo,q) one has 

IW ? ) = 2i J ^3 tr [M n G(p)T m (p,p + q)G(p + q)} , 

(12) 

where the mass vertex T m is a 4x4 matrix (the sublat- 
tice/valley index is omitted for clarity), and the factor 
of 2 accounts for spin. In the ladder approximation T m 
satisfies the self-consistent equation 



F m (p, q) = M m + ie' 



d 3 k G{k)T m (k,q)G{k + q) 



where (we set vf = 1 henceforth) 

fco + oik 



G(k) 



k 2 , — k 2 + ie 



(13) 



(14) 



To solve this set of equations, it is convenient to decom- 
pose r m in a basis of 4x4 matrices with well defined trans- 
formation properties under the SU(2) valley symmetry. 
Defining M = T z cr z , this basis may be taken as the four 



matrices M,T,a l which are scalars under this symme- 
try, and the matrices M ni T„, a l T n , which transform like 
a spin 1/2. With this choice we express T m as 



t m M + f ° I + TLa 1 



r m 



r° t 



(15) 



The equations are further simplified when T\ lm is ex- 
pressed in terms of its longitudinal and transverse parts 



q ■ r r 



T 1 



q x r r - 



(16) 



where q = q/q, and a similar relation applies for 
With the identities 



k ■ T nm — k ■ q T nm k x q T nn 



k x r, ; 



k-qTi 



kxqT* 



(17) 
(18) 



substituting Eq. (|15l) into Eq. (jT2"j) . and performing the 
trace, we obtain 



d 3 p 



T 
am 



[/nr nm + fuT 

/l4F° m )] , 



(19) 



(2tt) 3 D 
+px q(f 13 T^ 

where we have defined the denominator 

D(p,q) = [p 2 ~f + ie][(p + qof-(p + q) 2 +itl (20) 

and where all (p, q) (specified below) are even func- 
tions under the reversal of the relative angle 6p t $ = 
6p — q . Because of the decomposition in Eq. (|15[) . the 
scalar parts decouple completely and are not needed. We 
can then obtain equations for the remaining components 
of T m by multiplying Eq. (|13p by the corresponding basis 
matrices and taking the trace. One then obtains 



r„ 



Snm - ifi 



d 3 k 1 1 



(2tt) 2 £> £ 

+ k X q(fl3^n m - fli^nm) 

d 3 k 1 1 



[/nr. 



■n in + /i 2 r 

n i n. 



-0 



(2^) 2 £|p*_fc 

+k x q(f 23 r L nm + / 24 r° m )) 



[f2lT r , 



/ 22 r 



T 

n ni 



(21) 



(22) 



T mn and T mn satisfy similar equations, but arc not 
needed in what follows. We now perform a circular har- 
monic expansion 



p(« P ,n q ) 



dO 



2ir 



2tt 



r(p,g), (23) 



and retain only the first order contribution. Terms con- 
taining k x q are odd and vanish. Thus, and T° 
completely decouple to first order. Moreover, from the 
structure of Eqs. (jT9")) , (|21[) and (|22l) it can be seen that 



in fact n„ m = 8 nm Tl. As expected, the Kekule (IIn and 
II22) and CDW (II33) response functions are the same. 
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With this simplification the relevant components of /j. 
are 



fi 



-k (k + q ) + k(k + if), 



J 12 = /2l = H— k q), 

2(qx k) 2 -> -> 

f-22 = a + fc (fco + go) - + g). 

Defining the dimensionless kernels 



and 



C (n) (a;) 



d9. 



2ir (1 + x 2 + 2xcos6 k ) 1 / 2 ' 



(24) 
(25) 

(26) 
(27) 
(28) 



the self-consistent equations to first order in the circular 
harmonic expansion finally read 



r(0 ,o) = i + £ f dkcW(K[¥r(°>°) + KgrpV), 

2p J 



(29) 
(30) 



A numerical analysis shows that the mixing kernel K\i 
is small compared to Kn and may be neglected also. In 
this case, the final equations determining the response 
function, spelling all momentum dependence, read 

rW(p,q) = 1 + L [ A dkC^(k/p)K{f(k,q)T^(k,q), 
Z V Jo 

(31) 

U(q) = - f dpK^{p,q)T^){p,q), (32) 
~x Jo 

where A is an ultraviolet cutoff regularizing the inte- 
grals. Note the product {k/p)K^ \k,q) goes as p/k 
for large k, so the iteration of this equation produces a se- 
ries of logarithms characteristic of power law behaviour. 
Also note that when the external g < go, all develop 
an imaginary part for (qo — q)/2 < k < (go + g)/2. 

We solve Eq. (|3T|) numerically by discretizing the mo- 
mentum k on a logarithmic mesh and solving the corre- 
sponding matrix equation by Gaussian elimination. The 
integration of Eq. (|3"2"j) is straightforward. The result 
of this procedure is Il(go,g)- It is convenient to repre- 
sent it as the difference All = n(q ,go + Sq) — n(q ,q ) 
with Sq = go — g. Fig. [2] displays the real and imag- 
inary parts of All for \Sq\ « q . We observe a cusp 
at 6q = in the real part, and a finite imaginary part 
for 6q < 0. Log plots of both sides of the real part and 
and the imaginary part reveal power laws as Sq — > 0. A 
Kramcrs-Kronig analysis for \Sq\ << go shows that this 
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FIG. 2. Mass susceptibility AII(go,go + Sq) for \8q\ « qk 
and j3 = 0.36, real part (full line) and imaginary part (dashed 
line). 
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FIG. 3. Logarithmic plot of the mass susceptibility 
AH(q , qo + Sq) for P = 0.24 and Sq > (dotted line). The full 
lines are linear fits with r\ = 0.45 for Sq « qn and r/o = 0.69 
for Sq » qK- Inset: The exponent rj as a function of p. 



is only consistent if All cx (Sq) v , i.e. the exponents are 
all the same. Fig. [3] shows a log plot for 5q > where 
power law behavior is evident for Sq << qo. We also ob- 
serve that All crosses over to a different power law for 
Sq >> qo, which we identify as the static result q no — . The 
inset of Fig. [3] shows that rj is /3-dependent, and that it 
tends to the non- interacting result in Eq. (fTU)) as B — » 0. 
The dependence of 770 on /? can be found in Ref. 9. 

Finally, we plot the phonon dispersion relation, which 
is the main result of this work. This is given in terms 
of the self-energy evaluated at the phonon frequency ujk- 
To ease the comparison at different values of B, we will 
also represent the difference 



A 



Aw(g) = u(q) -u(q K ) = ^- (IL(u K ,q) - ^(^k,Qk)) 

(33) 

where we have recovered physical units with hvp = 6.5 
cVA The values of the parameters used are Xk =0.1 
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and A = 1.7eV. The phonon dispersion is depicted in 
Fig. @] for different values of j3. The dispersion follows 
the static power law q r > o( - 13 ^ for q » qjc, and the cusp 
turns into q v ^ as discussed. 




0.00 0.05 0.10 0.15 0.20 0.25 
q(A"') 

FIG. 4. Ai phonon dispersion relation Acj(q) measured from 
the K point for ft = 0,0.1,0.2,0.3, with higher curves cor- 
responding to higher values of /?. Note that ui(qK), which 
depends on /3, has been substracted from each curve for an 
easier comparison. Inset: the Kekule phonon displacements. 



IV. DISCUSSION 

Our computation has shown that interactions turn the 
Kohn anomaly at the K point into a power law, so it is 
natural to ask whether the same effect happens for the 
anomaly at T. This is not expected in general grounds, 
because the corresponding self-energy is built with ver- 
tices corresponding to a conserved current, and these 
type of operators do not have anomalous dimensions be- 
cause of Ward identities^. This is also consistent with 
the fact that the Coulomb interaction renormalizes the 
Ai electron-phonon coupling strongly, but not the E2 
on o 21 i 23 . Power law behavior is thus only expected in 
the K point anomaly. 

From the experimental point of view, there are sev- 
eral techniques available for the measurement of the A\ 
phonon dispersion, and each one has its own potential 
difficulties. In general, the power law at q > qx appears 
in a range of momenta that has been already probed with 
different techniques, while the cusp structure lies within 
the precision limits of current experiments, and may re- 
quire more effort. 



Electron Energy Loss Spectroscopy (EELS) is for 
example a suitable technique that has already been 
used to map the phonon dispersion at the K-point in 
graphene. This experiments have been performed on dif- 
ferent substrates for which graphene behaves as quasi- 
freest andingi^^, such as Pt (this is important as hy- 
bridization with the substrate strongly changes the elec- 
tron band structure and the Kohn anomaly^!). Metallic 
screening is however a disadvantage as it spoils the crit- 
ical behaviour of the electrons, and an insulating sub- 
strate would be more suited to observe the effect. 

A more indirect experiment (with insulating substrate) 
is to track the dependence of the 2D Raman peak with 
incoming laser energy. This method has been usecU^ to 
measure the dispersion of the A\ phonon. While the 
amount of data it yields and the range of momenta it cov- 
ers is limited and not very close to the K point, the obser- 
vation of the q > qk regime is certainly possible. Finally, 
X-rays are a usual tool to measure phonon dispersions 
in 3D crystals, and while it is probably challenging to 
obtain enough intensity from a single sheet of graphene, 
experiments in graphit o 15 ' 32 might be used to deduce the 
phonon dispersion. This approach is not straightforward 
because the electronic structure of graphite is different 
from graphene, and this must be taken into account. 
Nevertheless, it is encouraging to observe that precision 
measurements show an A\ phonon dispersion that is not 
at all linear—. 

A final comment concerns the robustness of our result 
to more refined approximations schemes than the ladder 
summation. While other sets of diagrams may modify 
our quantitative predictions, it is very unlikely that the 
non-analytic behaviour can be removed in this way. One 
may consider, for example, the inclusion of self-energy 
terms for the electron propagator™, which may produce a 
slow logarithmic dependence of the exponent. Finally, we 
also note that the 1/N approximation does gives power 
law behaviour for the Kekule mass correlato r 10 ' 11 (and 
thus the self-energy) as well. 

In summary, this work has shown that the elusive crit- 
ical behavior of interacting Dirac electrons in graphene 
manifests itself through a power law Kohn anomaly for 
the A\ phonon at the K point. 

V. ACKNOWLEDGMENTS 

We thank A. Politano for very useful discussions. 
Support from NSF through Grant No. DMR- 1005035, 
and from US- Israel Binational Science Foundation (BSF) 
through Grant No. 2008256 is acknowledged. 



1 A. H. Castro Neto, F. Guine a, N. M. R. Pere s, K. S. 
Novoselov, and A. K. Geim, |Rev. Mod. Phys.| 81, 109 
(2009) 



2 V. N. Kotov, B. Uchoa, V. M. Pereira, A. H. Castro Neto, 
and F. Guinea, Rev. Mod. Phys., submitted 

3 J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, 
|Nucl. Phys. B] 424, 595 (1994) 



6 



4 D. E. Sheehy and J. Schmalian, |Phys. Rev. Lett.| 99, 
226803 (2007) 

5 D. V. Khveshchenko, |Phys. Rev. Lett. 1 87, 246802 (2001) 

6 O. V. Gamayun, E. V. Gorbar, and V. P. Gusynin, 
|Phys. Rev. B] 81, 075429 (2010) 

7 J. P. Reed, B. Uchoa, Y. I. Joe, Y . Gan, D. Casa, E. Frad- 
kin, and P. Abbamonte. ISciencel 330. 805 (2010) 

8 D. C. Elias, R. V. Gorbachev, A. S. Mayorov, S. V. Mo- 
rozov, A. A. Zhukov, P. Blake, L. A. Ponomarenko, I. V. 
Grigorieva, K. S. Novoselov, F. Guinea, and A. K. Geim, 
Nat. Phys.(2011) 

9 J. Wang, H. A. Fertig, and G. Murthy, |Phys. Rev. Lett.| 
104, 186401 (2010) 

10 J. Gonzalez, [Phys. Rev. B| 82, 155404 (2010) 

11 A. Giuliani, V. Mastropietro, and M. Porta, |Phys. Rev. B| 
82, 121418 (2010) 

12 M. Franz, T. Pereg-Barnea, D. E. Sheehy, and 
Z. Tesanovic, [Phys. Rev. 5] 68, 024508 (2003) 

13 V. P. Gusynin, D. V. Khveshchenko, and M. Reenders, 
|Phys. Rev. B[ 67, 115201 (2003) 

14 D. L. Mafra, G. Samsonidze, L. M. Malard, D. C. Elias, 
J. C. Brant, F . Plentz, E. S. Alves, and M. A. Pimenta, 
|Phys. Rev. B[ 76, 233407 (2007) 

15 A. Griineis, J. Serrano, A. Bosak, M. Lazzeri, S. L. 
Molodtsov, L. Wirtz, C. Attaccalite, M. Krisch, A. Ru- 
bio, F. Mauri, and T. Pichler, |Phys. Rev. 13] 80, 085423 
(2009) 

16 A. Politano, A. R. Marino, V. Formoso, and G. Chiarello, 
ICarbonl 50. 734 (2012) 

17 A. Luther a nd I. Peschel,|Phys. Rev. B| 9, 2911 (1974) 

18 C. Chamon, |Phys. Rev B| 62, 2806 (2000) 

19 H. Suzuura and T. Ando, J. Phys. Soc. Jpn. 77, 044703 
(2008) 



20 A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, 
M. Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. 
Novoselov, S. Roth, and A. K. Geim, |Phys. Rev. Lettl. 97, 
187401 (2006) 

21 M. Lazzeri. C. Attaccalite. L. Wirtz. and F. Mauri. 
|Phys. Rev. 5] 78, 081406 (2008) 

22 S. Viola Kusminskiy, D. K. Campbell, and A. H. Cas- 
tro Neto, |Phys. Rev. B] 80, 0354 01 (2009) 

23 D. M. Basko and I. L. Aleiner, [Phys. Rev. B| 77, 041409 
(2008) 

24 D. M. Basko, S. Piscanec, and A. C. Ferrari, [Phys. Rev. B[ 
80, 165413 (2009) 

25 S. Piscanec, M. Lazzeri, F. Mauri, A. C. Ferrari, and 
J. Robertson, |Phys. Rev. Lett. 1 93, 185503 (20 04) 

26 M. Lazzeri and F. Mauri, |Phys. Rev. Lett.| 97, 266407 
(2006) 

27 A. H. Castro Neto and F. Guinea, |Phys. Rev. B| 75, 045404 
(2007) 

28 W.-K. Tse, B . Y.-K. Hu, and S. Das Sarma, 
|Phys. Rev. LetEl lOl, 066401 (2008) 

29 S. Pisana, M. Lazzeri, C. Casiraghi, K. S. Novoselov, A. K. 
Geim, A. C. Ferrari, and F. Mauri, INat. Materl 6, 198 
(2007) 

30 H. Yanagisawa, T. Tanaka, Y. Ishida, M. Mat- 
sue, E. Rokuta, S. Otani, and C. Oshima, 
|Surface and Interface Analysis 37, 133 (2005) 

31 A. Allard and L. Wirtz. INano Letters! 10, 4335 (2010) 

32 J. Maultzsch, S. Reich, C. Thomsen, H. Requardt, and 
P. Ordej6n, [Phys. Rev. Lettl. 92, 075501 (2004) 



